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ABSTRACT 


Two-dimensional  windows  find  applications  in  many  diverse  fields, 
such  as  the  spectral  estimation  of  random  fields,  the  design  of  two- 
dimensional  digital  filters,  optical  apodization,  and  antenna  array 
design.  Many  good  one-dimensional  windows  have  been  devised,  but 
relatively  few  two-dimensional  windows  have  been  investigated.  In  this 
paper  we  show  that  good  two-dimensional  windows  can  be  obtained  by 
rotating  good  one -dimensional  windows.  That  is,  if  w(x)  is  a  good 

r~ 2  2~ 

symmetrical  one-dimensional  window,  then  W2(x, y)  =  w(\x  +y  ) 
is  a  good  circularly  symmetrical  two-dimensional  window. 


Accepted  for  the  Air  Force 

Joseph  R.  Waterman,  Lt.  Col.,  USAF 

Chief,  Lincoln  Laboratory  Project  Office 


iii 


TWO-DIMENSIONAL  WINDOWS 


I.  Introduction 

Two-dimensional  windows  find  applications  in  many  diverse  fields,  such  as 
the  spectral  estimation  of  random  fields,  the  design  of  two-dimensional  digital  filters, 
optical  apodization,  and  antenna  array  design. 

1-4 

Many  good  one-dimensional  windows  have  been  devised,  however,  rela¬ 
tively  few  two-dimensional  windows  have  been  investigated.  ^In  this  paper,  we 
establish  a  result  which  enables  us  to  get  good  two-dimensional  windows  from  good 
one-dimensional  windows. 

II.  The  Problem 

We  first  review  briefly  the  one-dimensional  problem.  Let  the  Fourier  trans¬ 
form  of  a  function  f(x)  be  F(u).  For  some  reason,  we  want  to  truncate  f(x): 

g(x)  =  f(x)  w(x)  (1) 

w(x)  =  0,  for  | x  |  >  A  (2) 

where  A  is  a  constant. 

The  Fourier  transform  of  g(x)  is 

G(u)  =  F(u)<8>  W(u)  (3) 

where  W(u)  is  the  Fourier  transform  of  w(x),  and  ®  denotes  convolution. 

Our  problem  is  to  choose  an  appropriate  shape  for  the  window  function  w(x)  such 
that  G(u)  is  close  to  F(u)  and  in  any  region  surrounding  a  discontinuity  of  F(u), 

G(u)  will  not  contain  excessive  ripples.  It  is  well  known  that  the  Fourier  trans¬ 
form  W(u)  of  a  good  window  w(x)  should  have  a  big  central  peak  and  small  sidelobes. 
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In  two  dimensions,  the  problem  is  entirely  similar.  Let  the  two-dimensional 
Fourier  transform  of  a  function  f2(x,y)  be  F2(u,v),  and  let 


where 


g2(x>y)  =  f2(x>y)  w2(x,y) 
w  2(x,  y)  =  0  for  |  x2  +  y2 1  >  A2 


(4) 

(5) 


The  Fourier  transform  of  g2(x,y)  is 


G2(u,v)  =  F2(u,v)<g)  W2(u,v) 


(6) 


where  W2(u,v)  is  the  Fourier  transform  of  w2(x,y).  The  problem  is  to 
choose  an  appropriate  shape  for  the  two-dimensional  window  function  W2(x,y)  such 
that  G2(u,v)  is  close  to  F2(u,v)  and  in  the  neighborhood  of  a  discontinuity  of  F2(u,v), 
G2(u,v)  does  not  contain  excessive  ripples. 

III.  The  Result 


Intuitively,  we  feel  that  if  w  (x)  is  a  good  symmetrical  one -dimensional 

window,  then  , - 

w2vx,y)  =  w  (x/x+y2)  (7) 


will  be  a  good  two-dimensional  window. 

This  is  indeed  partially  verified  by  the  following  two  examples. 


The  first  example  is 


w  (x)  = 


1,  for  |  x  |  £  1 
0,  for  |  x  |  >  1 


(8) 


Then  the  Fourier  transform  is 

2  sinu 


Wfu)  = 


u 


(9) 


whose  first  side-lobe  peak  is  about  20#of  the  peak  at  u  =  0.  The  corresponding 
two-dimensional  window 
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!2  2 
1,  for  |  x  +  y  |  s  1 

0,  for  |  x^  +  y2  |  >1 

has  the  two-dimensional  Fourier  transform 


W>,v)  = 


2TTJJ  (Vn2  +  v2) 


fl  2 

v/u  +v 


whose  first  side-lobe  peak  is  only  about  124  of  the  peak  at  u  =  0  =  v. 
The  second  example  is 

w  (x)  =  j  1  -  x,  for  |  x  |  £  1 
(0,  for  |  x|  >  1 

The  Fourier  transform  is 


W(u)  = 


4) 


(10) 


(ii) 


(12) 


(13) 


whose  side -lobe  peak  is  about  \<4  of  the  peak  at  u  =  0.  The  corresponding 
two-dimensional  window 


9  9  9  9 

w2(x,y)  =  j  1  -y/x  +  y  ,  for  |  x  +  y  |  ^  1 
0,  for  |  x2  +  y2 1 >  1 


has  a  two-dimensional  Fourier  transform 

P 

W2(u,v)  =  2n  "3 


I  J 


o(t)dt  -  p  JQ(p) 


/  2  2 

where  p  =\/u  +  v  ,  whose  first  side-lobe  peak  is  only 


(14) 


(15) 

about  2# of  the  peak 


at  u  =  0  =  v. 
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The  above  comparisons  are,  however,  unfair.  Because  when  we  convolve 
a  window  with  a  discontinuity,  what  count  in  the  one-dimensional  case  are  the 
areas  under  the  side-lobes,  while  in  the  two-dimensional  case  what  count  are  the 
volumes  under  the  side -lobes.  A  fair  comparison  would  be  to  look  at  the  result 
of  the  convolution.  One  thing  we  can  say  along  this  line  is  contained  in  the  following 
theorem  which  is  the  main  result  of  this  paper. 

Theorem.  If  a  symmetrical  one-dimensional  window  w  (x)  and  a  two-dimensional 
window  W2(x,y)  are  related  by 


(16) 


then  their  Fourier  transforms  W(u)  and  W^u.v)  satisfy  the  relation 


■jn  w2(u,v)  ®  H2(u,v)  =  W(u)®H(u) 


(17) 


where 


(18) 


H2(u,v)  = 


1,  for  u  sO  and  all  v 
0,  for  u  <0  and  all  v 


(19) 


and®  denotes  convolution. 


Proof.  We  first  show  that 


CO 


(20) 


By  definition, 


CO 


wo(x,y)  =  — r- 
z  4rr 


j(xu+yv) 


CD 
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Whence 


w2(X 


du  e' 


jxu 


b  f  dv  W2(u,v) 


But  from  Eq.  (16), 

w2(x,0)  =  w(x) 

Therefore  from  Eq.  (21),  dv  W2(u,v) 

”00 

is  the  Fourier  transform  of  w(x).  This  established  Eq.  (20). 

Now,  00 

W(u)  <g>  H(u)  =  f  dt  W(t)  H  (u  -t) 


(21) 


and 


-I 


dt  W(t) 


W2(u,v)  ®  H2(u,v)  =  If  dtds  W(t,s)  H2(u-t,v-s) 


U  CO 

=  J  dtj  dsW(t.s) 


From  Eqs.  (22)  and  (20),  we  have 

U  03 

W(u)®  H(u)  =  fttf  ds  W(t,  s) 


“00  “00 


1 


=  2n  W2(u,v)®  H2(u,v) 
by  virtue  of  Eq.  (23). 


(22) 


(23) 
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IV.  The  Design  of  Two-Dimensional  Non-recursive  Ideal  Low-Pass  Filters 


One  way  of  designing  one -dimensional  non-recursive  digital  filters  is  the 
7 

so-called  window  method  .  To  fix  ideas,  let  us  consider  the  design  of  an  ideal 
low-pass  filter.  The  ideal  frequency  response  is 


F(u)  =  j  1,  for  |u|  £B  (24) 

|  0,  for  |u  |>  B 

where  B  is  a  constant.  (This  frequency  response  is  actually  repeated 
periodically  because  the  impulse  response  is  sampled).  The  inverse  Fourier 
transform  f(x)  of  F(u),  which  is  the  impulse  response,  has  infinite  duration,  but 
in  reality  we  have  to  use  a  finite-duration  impulse  response.  So  we  truncate  f(x) 
by  using  a  window: 

g(x)  =  f(x)  w(x)  (25) 

w(x)  =  0,  for  | x |>  A  (26) 

where  A  is  a  constant.  The  actual  frequency  response  we  are  getting  is 

then 

G(u)  =  F(u)®  W(u)  (27) 


Suppose  now  we  wish  to  design  a  two-dimensional  ideal  low-pass  filter  with 
an  ideal  Frequency  response 


F2(u,v)  =  F( 


+  v2)  = 


1,  for  |u2  +  v2 1  £  B2 
0,  for  |u2  +  v2|  >B2 


(28) 


And  we  truncate  the  two-dimensional  impulse  response  f2(x,y)  by  a 

dimensional  window  w  (x,y).  Then  the  actual  frequency  response  we  are 

-2 

G2(u,v)  =  F2  (u,v)<&  W2(u,v) 


two- 

getting  is 
(29) 
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Let  us  assume  that  the  widths  of  W(u)  and  W2(u,v)  are  much  smaller  than 
B,  then  near  the  discontinuities  of  F  (u),  viz. ,  u  =  ±  B,  G(u)  is  essentially  equal  to 
the  convolution  of  W(u)  and  a  one-dimensional  step  function,  and  similarly,  near 


2  2  2 

the  discontinuities  of  F2(u,v),  viz.  ,  u  +  v  =  B  ,  G2(u,v)  is  essentially  equal  to 


the  convolution  of  W^u.v)  and  a  two-dimensional  step  function.  It  therefore  follows 
from  our  theorem  that: 


if 


(30) 


then 


(31) 


This  means  that  we  can  design  a  good  two-dimensional  low-pass  filter  by 
using  the  window  w2(x,y)  as  given  by  Eq.  (30),  if  w(x)  is  a  good  window  to  use 
in  designing  a  good  one-dimensional  low-pass  filter. 

V.  Summary 

We  have  shown  that  if  w(x)  is  a  good  symmetrical  one-dimensional  window, 
/ 2  2 

then  W2(x,y)  =  w  (>/x  +  y  )  is  a  good  circularly  symmetrical  two-dimensional 
window. 
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